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1. Introduction 

Semiconductor quantum dots are often referred to as "artificial atoms" owing to their 
discrete energy spectrum. The progress in nanotechnology has made it possible to 
employ quantum dots as a sui generis solid state laboratory for studies of quantum 
mechanics pQ. Interband optical pumping of quantum dots gives rise to the electron- hole 
pairs or zero-dimensional excitons, which, as shown in the pioneering works [21 [3l H], 
can strongly couple with the photon trapped in the optical microcavity. The strong 
coupling effect results in a coherent energy transfer between the photon and the 
exciton, this phenomenon is widely studied for bulk materials and planar quantum 
well structures [5j El [7J [8]. Its observation in zero- dimensional systems has attracted 
an enormous excitement of the research community since the concepts of quantum 
electrodynamics were directly transferred to the solid state. Such quantum-dot-in-a- 
cavity systems demonstrate fascinating fundamental physics [91 [TUl [Til [12] and may be 
advantageous for quantum optics device applications [131 [HI [HI EE] • 

The physical concept of the strong coupling in zero- dimensional microcavities can 
be easily understood considering, for the sake of an example, two classical pendulums 
with close frequencies of oscillations U\ and uj<i connected by a spring, which induces the 
coupling between the pendulums of a strength g. One of these oscillators represents 
a photon, the other one stands for an exciton, and the spring describes radiative 
recombination of the exciton into the photon mode. If the dampings of the individual 
oscillators are small compared with the coupling constant, the eigenfrequencies u of this 
interacting system follow from the simple equation 



The normal modes correspond to the coupled oscillations: By exciting one pendulum 
one eventually excites another, so the energy is transferred back and forth between 
them. 

The analogy between the quantum electrodynamical problem of the quantum dot 
exciton interacting with the microcavity photon and purely classical problem of two 
coupled pendulums is quite deep. The excitonic polarization in a semiconductor is 
indeed described by the oscillator-like equation of motion and so does the electric 
field of the cavity photon [U [7J. Hence, the physics of these two different systems 
is similar which greatly simplifies theoretical description of the quantum-dot-in-a-cavity 
dynamics [161 The situation becomes particularly interesting if AT > 1 dots are 
placed in the microcavity: It turns out that, provided the dots are identical, only 
one excitonic mode - termed superradiant - interacts with the cavity mode, and the 
interaction constant is enhanced compared with that in one dot by the factor \/N. In 
classical language, for iV pendulums describing quantum dot excitons, one oscillation 
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mode is specific, namely, the mode where all pendulums oscillate in phase with each 
other. 

Just like the classical oscillators which are ideal only as far as model situations are 
considered, the excitons and photons interact with the environment, which gives rise to 
their damping and dephasing, moreover, the oscillation law can differ from the harmonic 
one. Main reasons of unharmonicity are identified: interactions between excitons lead 
to their energy shifts, the oscillator strength saturation results in the decrease of the 
coupling constant with an increase of the exciton number in the system, excitons can 
bind together to form biexcitons [HUH]. Hence, the analysis of the non-linear dynamics of 
the quantum dot excitons strongly coupled to the cavity mode is of prime importance. 
In particular, the stability of the superradiant mode with respect to the interactions 
should be investigated. 

In this paper, we consider the simplest possible and physically most transparent 
case of the non-linearity caused by the exciton-exciton interaction. It is analogous to the 
cubic unharmonicity of the pendulum^ We calculate the optical emission spectra of the 
quantum-dot-in-a-cavity system under a non-resonant excitation (photoluminescence 
spectra). The non- resonant pumping is modeled as a random force acting on the 
corresponding oscillator [T7J [22]. The parameter which controls the non-linearity is 
the pumping rate: the higher the pumping, the larger the fluctuations of excitonic 
polarization and, correspondingly, the higher are the non-linear terms in the equations 
of motion. 

In our approach, the quantum dot state is described as a classical unharmonic 
oscillator. Hence, this approach is valid for large enough quantum dots where excitons 
are quantized as a whole. Similar behavior can be expected in a variety of systems 
including the semiconductor systems with planar quantum microcavities where excitons 
can be trapped by the disorder and the non- linear regime is easily reached p3j [23] , 
or quantum well structures with dipolar excitons j25j EE], as well as many other, e.g. 
optomechanical structures where the optical mode of the cavity interacts with a classical 
oscillator [27]. 

The paper is organized as follows. Section [2] outlines the model, Sec. [3] is devoted 
to the role of the reservoir fluctuations, Sec. H] presents the study of the nonlinearity due 
to the exciton interaction. Main results are summarized in Conclusions, Sec. EJ 

2. Model 

For distinctness we consider here a zero-dimensional microcavity where one or several 
quantum dots are embedded. We suppose that the energy (or frequency, we put h = 1 
for brevity) spacing between the cavity modes is large enough to consider only one 
relevant photonic mode whose frequency uc is close enough to frequencies of optical 
transitions in quantum dots ux,i, i = 1, ...,iV, where N is the dot number. For 

| The effects of the oscillator strength saturation related with the two-level nature of exciton transition 
in quantum dot were considered in Refs. [T§1 [201 \2T\ . 
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simplicity, the interaction of the cavity mode with the ground states of quantum dot 
excitons is considered only, the treatment can be generalized to allow for the excited 
states. Moreover, the polarization degree of freedom of the cavity mode and spin degrees 
of freedom of excitons are disregarded. Under these approximations the equations of 
motion for the dimensionless electric field E and excitonic polarizations Pi can be written 
as 

dE 



dt 

dp f. r 



(3) 



(// Au x> i + ^ + iai\Pi\ 2 + ipin R )Pi-igE + fi(t), i = l...N. 

Here Tc and r^, are (phenomenological) decay rates for the cavity and excitons, 
respectively, g is the light-matter coupling constant (taken the same in all dots for 
the sake of simplicity). Its evaluation is beyond the current work, rigorously it can be 
done by solving Maxwell equations and Schrodinger equation for the exciton envelope 
functions in quantum dots, see Refs. [281 [29J [30] . 

Equation (jHJ) for exciton polarization contains also non-linear and driving terms. 
Former ones, proportional to the interaction parameters > 0, describe the blueshifts 
of excitonic states due to interactions within the same dot. Here and further we 
assume that the quantum dot size is comparable or larger than excitonic Bohr radius 
to accomodate several excitons in the dot |19jJ§l Terms cx ^Ur take into account 
the interaction of quantum dot excitons with a reservoir formed, e.g. by excitons 
and electron-hole pairs within the wetting layer. The reservoir population ur is, 
in general, function of the pumping intensity W and exciton occupations, ur = 
ur(W, \Pi\ 2 , ■ ■ ■ , l-P/vl 2 )- Driving terms described by random forces f, (f) = account 
for the exciton generation in quantum dots caused by their relaxation from the wetting 
layer and excited states, see Refs. (22J, [T7j for details. These random forces represent 
the white noise 

(fi(t)f*(t')) = SAAt-t'), (4) 

characterized by the exciton generation rate Si proportional to the reservoir population, 
Si = SiUR, where is related with the relaxation rate of excitons towards the ground 
state [17]. In what follows we assume that higher-order correlators of random forces are 
reduced to the second order ones, Eq. in accordance with the Gaussian distribution. 
Such random forces model the incoherent non-resonant pumping of our system. Similar 
approaches were used to study optical emission of exciton-polaritons in Bragg multiple 
quantum well structures and for planar microcavities [22j [32j [33], [34] . 

We are ultimately interested in the luminescence spectrum given (up to the common 
factor) by [T7] : 

/oo 
(E(t + t')E*(t))e^'dt', (5) 
-oo 

§ Another possibility is to consider a planar microcavity structure with a lateral potential confining 
the exciton, see e.g. Ref. [51] . 
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where averaging over time t is assumed. If the non-linear contributions are disregarded, 
emission spectrum can be obtained analytically with the result [17] : 

Here, fl m are the eigenfrequencies of the homogeneous system (jHJ) found neglecting non- 
linear terms oc a i} (3 i.e. Vt m are the polariton frequencies, and [Ci$ ,C\ , ...,C^\ = 
[E, Pi, . . . , Pn] are the corresponding eigenvectors, i.e., the Hopfield coefficients of the 
excitonic polaritons formed by the excitons coupled to the cavity mode [7J . The spectrum 
Eq. ([6]) can also be recast as a sum of terms with poles at polariton frequencies Q m , Q^. 

The detailed analysis of Eq. (jHJ) is presented in Ref. [IT] . Here we briefly discuss an 
important liming case where all exciton frequencies are the same, ux,i = wx,2 = • • • = 
ux,n = wi, exciton decay rates are the same, Txa = Tx,2 = • • • = ^x,n = Tx and the 
decay rates are negligible compared with the coupling constants, T X ,T C <C g- Under 
these assumptions the polariton frequencies and Hopfield coefficients take simple form. 
There are two mixed modes with the frequencies, cf. Eq. (j2J), 



n^ = — 1 —^ ± \J{-^) +N9 > M 

and all the remaining N — 1 modes correspond to the exciton states decoupled from 
light, fl m = ujx for m = 3, . . . N + 1. The effective coupling strength is enhanced by 
the factor y/~N due to the superradiance effect: dipole moments of excitons oscillate 
in phase. Correspondingly, emission spectrum has two peaks at Qi and fl 2 splitted by 

1 

(u - cox + iT x /2){uj- cue + iT c /2) - Ng 2 \ 2 
We note that Eq. (jSJ) is valid even if the condition Tx <C g or Tc g is violated [IT] . 
Two distinct peaks are observed in the emission spectra for the broadenings as high as 
Tx, Tc ~ y/Ng, otherwise these two peaks merge into one. 

Now we turn to the discussion of the non-linear effects. Firstly, we consider 
interactions with the reservoir described by the terms oc Pin^Pi, afterwards we discuss 
the effects of exciton-exciton interaction within the same dot described by the terms 

OC Ota I Pi I Pn . 



HnM = (lE^f) OC „ - /oW ,/,. , ^ /oX Ar „ 9|9 • (8) 



3. Interaction with reservoir 



Under conventional non-resonant pumping conditions the excitons are generated in the 
wetting layer and form a reservoir. If the pumping rate is moderate the majority of 
excitons are in the reservoir and its occupation is only weakly affected by the presence 
of the quantum dots, hence, one can treat ur in Eqs. (jBJ) as an independent quantity. For 
the steady-state pumping (n R ) = n R and there are certain fluctuations of 5n R around 
this time-averaged value, 8n R = n R — n R . 
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To elucidate the role of the reservoir we neglect completely the non-linearities caused 
by the exciton-exciton interaction in quantum dots, i.e. we put «j = in Eqs. fl3]). It is 
instructive to analyze the case of a single quantum dot and disregard the light-matter 
interaction. In this situation we obtain 

^ = - U)x + I y + ffiR + i/Wn«) P + /(*). (9) 

For clarity we have separated in Eq. (Q the contributions due to mean number of 
particles in reservoir and due to its fluctuations. Equation (jHJ) is the first order linear 
differential equation which can be integrated with the result 

P(t) = [ dt'f(t')ex P \-(iu x + i(3n R + r x /2)(t-t')-if3 [ 5n R (t")dt" 
Jo L Jf 

where the solution of homogeneous equation is omitted. The autocorrelation function 

(P(t)P*(0)), whose Fourier transform determines the single dot emission spectrum in 

the case of the regime of the weak coupling with the photon, reads 

1x .(io) 



(P(t)P*(0)) = A ( exp 



-iu x t ^ ifin R t -i/3 5n K (ii)dti 

* Jo 



Here the averaging over the random source realizations and reservoir population is 
assumed. The fluctuations of the reservoir particle number n R take place on much 
longer timescale compared with the fluctuations of random forces, Eq. (j3J). Hence, the 
averaging over the realizations of / is carried out independently. 

The oscillation spectrum, (|P(u;)| 2 ) = (P(t)P*(0))dt, which corresponds to the 
quantum dot emission spectrum in the weak coupling regime, depends strongly on the 
relation between various timescales in the system: exciton lifetime 1/Tx, correlation 
time in the reservoir r c , defined by the relation (5n R (t)5n R (0)) = (5n R )e~ t l Tc , and 
dephasing time caused by the reservoir, where 
1 f°° 

- = (3 2 dt (5n R (t)5n R (0)) = f3 2 (Sn 2 R )r c . (11) 

Td Jo 

If reservoir fluctuations are fast and small enough, r c r^, T^ 1 , the so-called motional 
narrowing regime can be realized and 

<^>l> K (,-^-^l[r, v /2 + 1 /(2 Tj) r < 12 > 

In this case interactions with reservoir slightly broaden the quantum dot spectrum. 

In the opposite case, where the fluctuations are strong and slow, the quantum dot 
spectrum can be presented as 

(|PH| 2 ) oc J dn R p(n R )5(u; - u x - Pn R ) = ip (j^j^) , (13) 

where p(n R ) is the distribution function of the reservoir, i.e. the probability that the 
number of particles in the reservoir is n R . The quantum dot spectrum Eq. f)13p in the 
regime of strong fluctuations is strongly asymmetric and reflects the distribution of the 
particles in the reservoir. Equation (IT51) is valid for the frequencies \uj — ux — 0n R \ > T x , 




Figure 1. (Color online) Effect of interaction with the reservoir (a) and within the 
quantum dot (b) on its spectrum (|P(w)| 2 ). Black/dashed, red/solid and blue/dotted 
curves in panel (a) correspond to linear regime, fast reservoir fluctuations regime and 
regime of slow reservoir fluctuations with Gaussian statistics, respectively. Calculation 
was performed at (5nR = 4Tx, l/r^ = 2Tx and j38fiR = 1.5Fx- Black/dashed, 
red/solid and red/dotted curves in panel (b) correspond to linear regime, exact 
solution, Eq. (|23|) . and approximate solution Eq. (|22|) . with (5-function replaced by 
Lorentzian with full-width at half-maximum equal to Tx ■ Calculation was performed 
at aS = 3Tx- 

that is why the spectral function of the quantum dot is replaced by 5-function in Eq. ( !T3|) . 
This equation also holds provided that the particle number fluctuations in the reservoir 
are so large that 

P^[^A)^>T X) t-\ (14) 

In this case the dephasing takes place at a timescale ~ (/3a/ (<5n^)) _1 and it is not 
sensitive neither to the oscillator lifetime nor to the correlation time of the reservoir. 

These two limiting cases of reservoir fluctuations described by Eqs. f|T2|) and ffT3l) 
are illustrated in Fig. [TJ^a) where the calculated quantum dot spectra are shown. 
Black/dashed curve presents spectrum for vanishing reservoir fluctuations, this spectrum 
is described by the Lorentzian with the full width at half-maximum (FWHM) equal to 
Tx, centered at the resonance frequency ux- Blue/dotted curve corresponds to the 
regime of fast reservoir fluctuations, Eq. (fl~2l) . This spectrum is also a Lorentzian, 
shifted to larger energies and broadened as compared to the linear regime. Red/solid 
curve is calculated in the regime of slow reservoir fluctuations. We have assumed, that 
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the reservoir population is characterized by a Gaussian distribution with given mean 
value fiR and dispersion Sur. In numerical calculations the ^-function in Eq. f lT3|) was 
replaced by the Lorentzian with FWHM equal to Tx- Then the convolution ( |T3|) yields 
Voigt distribution with Gaussian-like central part and Lorentzian wings, see Fig. Ufa). 

It is noteworthy that for the quantum dot placed into the microcavity, the emission 
spectrum in the limit of weak and fast fluctuations, r c <C r d , T x x is given by Eq. (jSJ) 
where ux is replaced by u x + [3n R and T x is replaced by T x + . The qualitative 
shape of the spectra in the case of strong or slow reservoir fluctuations can be obtained 
by the convolution of the reservoir distribution function p(nn) with Eq. ([6]) where ux is 
substituted by ux +f3riR. The resulting spectra are similar to those discussed in Sec. 14.21 
where the case of the nonlinearity within the dot is addressed. 

4. Interaction within the dot 

Now we focus on the non-linear effects due to the interaction of the excitons within the 
same quantum dot. Hereinafter we disregard fluctuations of the particles in the reservoir 
studied above in Sec. [3j Below we present one after another the studies of (i) the single 
dot case (Sec. 14.11) . (ii) the single dot coupled with the cavity mode (Sec. 14.2ft . and (iii) 
two quantum dots coupled with the cavity (Sec. 14.31) . 

4-1. Single dot 

The physical picture of the interaction effects within the dot on the emission spectra 
can be most transparently presented for the case of the single quantum dot which does 
not interact with the photon. Similarly to the situation studied in Sec. [3] its spectral 
function (\P(u)\ 2 ) defines the emission spectrum in the weak coupling regime. 
We start with the equation describing quantum dot as a non-linear oscillator: 



Similarly to the case of reservoir fluctuations, Eq. (Q, the interaction term in Eq. ffl5|) 
leads to the blueshift and broadening of the oscillation spectrum. Consequently, the 
analysis of reservoir fluctuations, performed above, may be used here. The strength of 
the fluctuating non-linear term is determined by the value of pumping. Large pumping 
(a(|P| 2 ) 3> Tx) corresponds to the regime of strong fluctuations, cf. Eq. ( fl4l) . Hence, 
the quantum dot spectrum may be presented in the form, similar to Eq. (1131) . 



where the distribution function F(Q) with the oscillator frequency Q = u>x + ct|-P| 2 , is 
determined from the statistics of the non-linear term in Eq. ( FTSj) . 

The major difference of the non-linear Eq. (115j) and the linear one Eq. describing 
the interaction of the quantum dot exciton with the reservoir is that the fluctuations of 
the polarization P themselves govern the blueshift and, on the other hand, the blueshift 




(15) 




(16) 
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determines the fluctuations of P. In order to determine the distribution function F(Q) 
we use the stochastic linearization of Eq. f]T5]) described in detail in Refs. [35] 136] . The 
starting point of the stochastic linearization is the time-dependent Eq. ( fl5l) where decay 
and pumping terms are neglected, T x = 0, f(t) = 0. The solution of Eq. (|T5|) is then 
given by P(t) = \P\e~ im , where 

n\P\ 2 = u x \P\ 2 + a\P\ A . (17) 

In the presence of the decay and pumping Eq. (JTTJ) does not hold. In the stochastic 
linearization approach the difference of the left- and right- hand sides of Eq. ( 1T7|) should 
be minimized. It implies, in particular, that the average (over the random sources 
realizations) blueshift of the resonance frequency is given by 

a(|P| 4 ) 

< n " = jp^y • ( 18 ) 

Hereinafter we use the notation (\P\ 2 ) = (\P(t = 0)| 2 ) = / doj(\P(oj) 2 \)/(27r). 
Equation f]T8]) is automatically satisfied when the distribution of the frequency Q is 
chosen in the form 

F(Q) = AT(Q - oj x )p (tt^L) , (19) 



a 

where M is the normalization constant, 

p(|p|2) = (^) exp ("(^y) ' (20) 

is the distribution function of the exciton polarization absolute value, and 

(I^IHt^, (21) 
J- x 

is given by the ratio of the pumping and decay rates. The shape of this distribution is 
independent of the pumping rate. In other words, the non-linear term in Eq. f|T5|) 
influences the oscillator phase only and does not affect its amplitude. Thus, the 
distribution Eq. ( 120]) is the same as in the linear regime and it inherits the Gaussian 
statistics of the noise term fluctuations. Using Eq. ([18]) and Eq. (120]) we find that 

. (u — OJx) exp — , u > ux 

(\P(oj)\ 2 ) oc F(oj) = 1 \ a(\P 2 \))' (22) 

0, oj < ojx- 

This function decays exponentially for large values of Q because the realization of 
high blueshift a|P 2 | is unlikely. Moreover, frequencies oj < ojx are not possible since 
interactions are repulsive and lead to the increase of energy only. Equation ( [22] ) 
demonstrates, that in the regime of high pumping the oscillator spectrum is strongly 
asymmetric and broadened due to the fluctuations of the resonance frequency. 

For arbitrary value of pumping strength the general analytical result for 
the oscillator spectrum can be obtained by means of the Fokker-Planck equation 
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technique [371 1ZE], similarly to the case of the noise-driven Duffing oscillator [35]. The 
spectrum reads 



In the linear in the pumping regime where S — > only the first term in the series 
Eq. (123"]) does not vanish and the result is reduced to Lorentzian with FHWM equal to 
Tx- For very large pumping, a(|P 2 |) 3> Tx, the series reduce to Eq. (1221) . 

The spectra of the single oscillator, (|P(o;)| 2 ), calculated for different pumping 
strengths, are shown in Fig. H^b). The spectrum at the strong pumping represented by 
red/solid curve, is strongly broadened and shifted as compared to the spectrum found 
in the linear regime. The latter is shown by black/dashed black curve. The width of 
the non-linear spectrum is of the same order as the blueshift. The spectrum found by 
the stochastic linearization, Eq. ( 122|) . and shown by red/dotted curve in Fig. DJb) well 
approximates the exactly calculated one, Eq. ( l23l) . The strong asymmetry of emission 
spectra calculated for the non-linear regime is clearly seen from the Figure. 

This concludes the discussion of the single oscillator and now we proceed to the 
analysis of the quantum dot coupled with the cavity mode. 

4-2. Single dot coupled with the cavity 

To start, we recall that in the linear-in pumping regime, S 0, and under the conditions 
of the strong coupling, g ^> Tc,Tx, the spectrum described by Eq. ([H]) with N = 1 
consists of the two distinct peaks at the frequencies of the system eigenmodes, excitonic 
polaritons. The main question we address here is how this two-peak spectrum changes 
with allowance for the non-linearity. In the general case the emission spectrum under 
strong pumping should be calculated numerically. This can be done either reducing the 
problem to the Fokker-Planck equation [37] or directly integrating the set of Eqs. (j3]). 
The latter procedure turns out to be more efficient, because the Fokker-Planck equation 
is computationally demanding already for N — 1 due to the large number of independent 
variables. In our calculations we have used the simplest 1.5-order Heun method for 
integration of stochastic differential equations [391 HO] • 

Luminescence spectra calculated for different detunings between the exciton and 
photon modes are presented in Fig. |2J Panel (a) shows the color plot of the emission 
intensity. Panels (b), (c), and (d) present the spectra for the detuning uc — ux equal 
to 6g, 0.5g and — 5g, respectively. Thin/red and thick/black curves are calculated, 




(23) 



where 
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Figure 2. (Color online) (a) Photoluminescence spectra for single quantum dot 
coupled to microcavity mode as function of the detuning between cavity and exciton 
frequencies. Panels (b),(c),(d) show the spectra for u>c = + Qg, = + O.bg, 
luc = <-0x — 5g, respectively. Black/thick and red/thin lines in panels (b), (c), and 
(d) correspond to non-linear and linear in pumping regime, aS = Q.3gTc and S = 0, 
respectively. These curves are normalized to their maximum values. Black/dashed 
curves in panel (a) are guides for eye demonstrating anticrossing. Calculations are 
carried out at Tx — O.lg, Tc = g, aS — 0.3gr<7. The time-domain integration was 
performed up to t = 100/g, with the time step 0.003/ g, the spectra were averaged over 
ensemble with 4000 realizations of random forces. 

respectively, for (i) the linear regime S — > and (ii) the regime where the non-linearity 
is already strong, S = 0.3Tc/<y. The spectra are normalized to their maximum values, 
the other parameters of the calculations are given in the caption to the Figure. The 
calculation demonstrates, that the non-linear spectra retain the characteristic two- 
peak structure, although the spectral shape is strongly affected by the non-linearity, 
in particular, it becomes asymmetric, as well seen in Figs. [2^b) and (d). The spectral 
maxima in the non- linear regime clearly exhibit the anticrossing behavior, see Fig. |2|a). 
From this we conclude that the strong coupling survives even for non-linear regime if 
*(\P 2 \)<9- 

Let us now discuss the spectra in Fig. [2] in more detail. We start from the case of 
large detuning between exciton and photon modes, \u x —^c\ ^ 9- in this situation the 
spectrum has two peaks related with the cavity and exciton emission. It is noteworthy, 
that in the non-linear regime, the blueshift of the exciton is determined by the dot 
population which, in its turn, is related with the detuning: (|P| 2 ) ~ S(ux— wc)V (Teg 2 ) 
(provided that Tx <C Tc and Tc <C \u>x — ojc\)- Hence, if bare exciton frequency ux is 
fixed, but the cavity frequency uc is varied, the position of exciton peak changes due 
to the variation of the blueshift. 

To confirm this argument we have plotted in Fig. [3] the dependence of the stationary 
intensities (\E\ 2 ) and (\P 2 \) on the detuning. For large detuning one has (\P 2 \) ^> (\E\ 2 ) 
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Figure 3. (Color online) Intensities of photon and exciton modes as functions of the 
detuning luc — ^x- Red/solid and black/dashed curves correspond to a(\E\ 2 ) /Tc and 
a(\P\ 2 ) /Tc, respectively. The parameters of calculations are the same as in Fig. [2] 



which is explained by the longer exciton lifetime as compared to the photon lifetime, 
Tc 3> Tx- For small values of detuning the curves become closer to each other as a 
result of the coupling between exciton and photon modes. 

Another important feature revealed in Fig. [2] is the strong asymmetry between the 
spectral shapes of the exciton peak for large positive and negative detuning, cf. Fig. [2(b) 
and Fig. [2(c). This is quite different from the linear regime, where the spectra Eq. (JBJ) 
are symmetric with respect to zero detuning. In the non-linear case the exciton peak 
is asymmetrically broadened due to the frequency fluctuations, like in the case of single 
oscillator, and the shape of the broadened peak depends on the sign of detuning. For 
negative detuning, uj c — uj x < 0, the high-energy tail dominates the spectrum similarly 
to the case of the dot decoupled from the cavity, see Sec. 14.11 and Fig. [1(b), while for 
positive detuning, uc — > 0, this tail is quenched. This is related with the fact 
that strong fluctuations of exciton polarization P are strongly suppressed since they 
correspond to large blueshifts where exciton mode approaches the cavity mode and, 
hence, short lifetimes. 

To understand this effect in depth we have plotted in Fig. |U the emission spectra 
[(a), (b)] along with the distribution functions of the quantum dot exciton intensity 
\P\ 2 . Distribution functions p(\P\ 2 ), shown by black/solid curves in Fig. 0(c) and 
Fig.H(d), were extracted from the numerical solutions of the system (JHJ). For comparison, 
red/dashed curves show exponential distributions (I20J) with the same average values 
of \P\ 2 - The emission spectrum of the microcavity can be presented in the following 
phenomenological form similar to Eq. (JT3J) and Eq. (fl~6j) obtained in the stochastic 
linearization method: 

J(w) oc J dttF(tt)I lin (cu) . (24) 
where In n (u) is given by Eq. (JHJ) with ux = ^ and the distributions F(Q) and p(\P 2 \) 




Figure 4. (a), (b) Cavity emission spectra calculated for the detuning wc — ^>x = 6.g 
(a), and ooc ~ u x = —5g (b). (c), (d) Distribution functions of the intensity \P\ 2 
calculated for the same parameters as in panels (a) and (b), respectively. Black/solid 
lines represent results of direct numerical calculation, red/dashed lines are calculated 
according to Eqs. ()20[) .(f24 |) . Other parameters of calculations are the same as in Fig. [2] 



are related by Eq. ([18]) . Corresponding spectra, calculated with numerically found 
functions p(\P\ 2 ), are shown in Fig. 11(a) and Fig. 11(b) by red/thin curves. We see, 
that Eq. f )24|) satisfactory describes the shape of the exciton peaks for both signs of 
the detuning and clearly demonstrates the correspondence between the shape of the 
exciton peak and the distribution function of the intensity \P\ 2 - For negative detuning 
the distribution function and the excitonic peak have exponential tails, similarly to the 
case of single oscillator, Eq. ( 1221) . Suppression of this tail for positive detuning means 
that the distribution function decays faster than exponential [cf. red and black curves 
in Fig. 11(c)] which is explained by interaction of excitonic and photonic modes. Indeed, 
for uc > ojx the non-linear blueshift of the exciton frequency decreases the detuning. 
Consequently, the exciton lifetime becomes effectively smaller due to the Purcell effect 
which suppresses the probability of such large detuning. Another effect leading to the 
same result is the repulsion of the excitonic mode with high blueshift from the cavity 
mode. In the opposite case, uc < ujx, the absolute value of the detuning is further 
increased by the blueshift, so the shape of the excitonic peak is not modified by the 
interaction with the cavity. 

Phenomenological Eq. (pMj) fails, however, to reproduce the ratio between the 
magnitudes of the excitonic and photonic peaks in the spectra. Formally, Eq. ( |24l) 
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Figure 5. (a) Photoluminescence spectra for two quantum dots coupled to 
microcavity mode as function of the detuning between cavity and exciton dot. Panels 
(b), (c), (d) show the spectra for u>c = ojx.i + 6g, cue = + 9, = — 6g, 
respectively. Black/thick and red/thin lines in panels (b) - (d) correspond to non-linear 
and linear in pumping regime, aS — 0.4grc and S — > 0, respectively. These curves 
are normalized at their maximum values. Black/dashed curves in panel (a) are guides 
for eye demonstrating anticrossing. Calculation was carried out for u)x,i — ^x.2 = .9, 
other parameters are the same as in Fig. [2l 



is valid provided that the timescale of the fluctuations of P is large compared with 
other timescales in the system [cf. Eq. ( [14]) ] or the magnitude of the blueshift exceeds 
by far all other energy scales. Neither of these conditions holds in the system under 
study. 

4-3. Two dots coupled with the cavity 

Now we turn to the discussion of the emission spectra of the cavity with two 
embedded dots, shown in Fig. |5j The specific feature of this problem in the linear 
regime is the formation of the collective, superradiant mode of the quantum dot 
excitons [HI [T7J [2TJ. In particular, as discussed above in Sec. El there are three 
eigenmodes of the homogeneous linear system Eq. (JBJ) for iV = 2. In case where 
&x,i — ux,2, one of these modes, with the excitons oscillating with the opposite phases, 
Pi = — Piy does not interact with the cavity mode, it is called dark. Two remaining 
modes correspond to the excitons, oscillating in phase, P\ = P2, and are formed by the 
coupling between the superradiant mode of the excitons and the cavity mode. Only 
two peaks are manifested in the emission spectrum. The superradiant effect leads 
to the enhancement of the Rabi splitting between these peaks from g to vNg, see 
Eq. (jSJ). When the frequencies of the dots ux,i are different, or the exciton tunneling 
between the dots is introduced, the dark modes mix with the superradiant mode and 
the spectrum acquires three-peak shape. As is demonstrated in previous work jTTJ , the 
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superradiant mode is stable against the disorder when the characteristic splitting of the 
dot frequencies is less than gy/N . Here we analyze the stability of the superradiant 
mode against the interactions. 

Figure |5] shows the emission spectra of the microcavity containing two quantum 
dots for various frequencies of the cavity mode. The frequency separation between the 
dots was fixed to be ux,2 -wi,i = g and the uc was varied. Overall intensity dependence 
on the emission and cavity frequency shown in Fig. [5(a) is quite similar to the intensity 
distribution for one dot in the microcavity [see Fig. 12(a)] and clearly demonstrates two 
peaks for any given cavity mode position. Thus, the superradiant mode is stable against 
the interaction. 

Panels (b), (c), and (d) of Fig. [5] present the details of emission spectra for three 
different cavity mode positions. Black/solid curves are calculated for the strong pumping 
where the nonlinear behavior is pronounced, red/thin curves correspond to the linear 
regime. The emission spectra in the non-linear regime corresponding to rather large 
detunings [panels (b), (d)] demonstrate the asymmetry studied for the single dot in the 
cavity in Sec. 14.21 Interestingly, the spectra in the linear regime demonstrate three peaks: 
two stronger ones correspond to the dot emission and weaker one to the cavity emission. 
The presence of non-linearity induced by the strong pumping qualitatively changes the 
emission spectrum: excitonic peaks merge and their intensity drops. Note, however, that 
the strong coupling regime is still maintained, see anticrossing in Fig. |5j The smallest 
splitting between the maxima of the spectra is larger than the Rabi splitting for the 
single dot 2g. This is a fingerprint of the superradiant mode stability in the non-linear 
regime. 

Thus, the superradiant mode becomes stabilized by interactions. Indeed, for 
positive detunings, uc > ^x,i,^x,2, the blueshifts of the quantum dots are different: 
The dot with higher exciton frequency (dot 2 in our calculation) has smaller lifetime 
due to proximity to the cavity mode, correspondingly, smaller exciton intensity (I-P2I 2 ) < 
(|Pi| 2 ), and, hence, smaller blueshift. Hence, with an increase of the pumping rate the 
blueshifted dot frequencies approach to each other, resulting in the decrease of the 
splitting between exciton frequencies and stabilization of the superradiant mode. 

5. Conclusions 

To conclude, we have developed a theory of non-linear emission of quantum dots coupled 
to the optical mode of the microcavity under non-resonant excitation. We model 
quantum dot excitons as non-linear oscillators taking into account the repulsive exciton- 
exciton interactions both within the dot and between the quantum dot and excitonic 
reservoir. We use random sources approach to model the relaxation processes and 
apply stochastic linearization and numerical integration of the Langevin equations to 
determine the spectra. 

Our model clearly shows that interactions: (i) blueshift the transition energy and 
(ii) to the same extent broaden the emission peak. The interactions result in the complex 
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behavior of the exciton lifetimes and intensities of emission as functions of pumping 
rate. The emission lines are strongly asymmetric and bear information on the exciton 
statistics. Contrary to the linear regime, the lineshapes are sensitive to the sign of the 
detuning between the exciton and photon modes. Interestingly, even for substantial 
pumping the strong coupling regime between the cavity mode and the quantum dot 
exciton can be preserved. Moreover, if two quantum dots are placed in the cavity, the 
superradiant behavior can be stabilized by the pumping. 
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